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EXECUTIVE  SUMMARY 


OBJECTIVE 

Develop  a  simple  model  for  HF  (2-30  MHz)  modem  performance  as  a  function  of 
various  channel  parameters.  Develop  techniques  for  predicting  those  parameters. 

RESULTS 

Modem  performance  was  modeled  based  on  three  parameters— signal-to-noise  ratio 
(SNR),  delay  spread,  and  Doppler  spread.  Each  of  these  parameters  is  measurable 
using  known  techniques.  The  model  demonstrates  the  inadequacy  of  using  only  SNR 
and  bit-error  rate  estimates  for  frequency  selection  in  automatic  link  establishment 
(ALE)  algorithms. 

RECOMMENDATIONS 

The  benefits  from  incorporating  delay  and  Doppler  spreads  in  the  ALE  frequency 
selection  algorithm  are  sufficiently  clear  to  justify  extension  of  the  algorithms  to 
include  these  factors.  Further  research  is  required  to  verify  the  sufficiency  of  SNR, 
delay  spread,  and  Doppler  spread  for  frequency  selection. 


Aooesslcn  For 

NXIS 

GRA4I 

DTIC 

TAB 

□ 

,  Unanin’inced 

n 

Juatl 

1  Icrt  !.cn _ 

Ry  .  ■ 

Diptrlbutior/ 

Avatloblltty  C’CiSua 

'»>r 

Dlst 

A'' 

I 

r 

1 

1 


CONTENTS 


INTRODUCTION .  1 

BACKGROUND  .  1 

OBJECTIVE  .  1 

APPROACH  .  1 

ANALYTICAL  ASSESSMENT .  2 

GAUSSIAN  NOISE .  2 

GAUSSIAN  NOISE  WITH  FLAT  FADING .  3 

EXTENSIONS  FOR  DQPSK .  5 

EXTENSIONS  FOR  FSK .  7 

ANALYTIC  ASSESSMENT  CONCLUSIONS .  12 

REFERENCES  .  15 

FIGURES 

1.  Normalized  rms-multipath  spread  approximation .  9 

2.  Bit-error  probability  for  16-tone  DQPSK .  13 

3.  Bit-error  probability  for  8-ary  FSK .  13 

TABLES 

1.  ALE  frequencies  and  bit  representations  .  10 

2.  Possible  error  occurrences  .  11 


INTRODUCTION 


BACKGROUND 

Standards  for  high  frequency  radio  automatic  link  establishment  (ALE)  are  being 
developed  at  the  federal  government  level,  within  the  Federal  Standards  series.  These 
same  standards  are  being  accepted  and  promulgated,  in  slightly  different  format,  by 
the  Department  of  Defense  as  NHL-STD- 188-141  A.  The  standards  support  automated 
means  for  selecting  the  “best”  operating  frequency  for  HF  communications  between 
two  locations.  The  standards  define  both  the  transmission  waveforms  and  frequency 
selection  algorithms. 

The  waveform  selected  for  the  ALE  standard  (8-ary  FSK)  is  not  the  same  as  either 
of  the  contemporary  waveform  types  used  by  the  Navy  (multitone  DQPSK;  equalized, 
high  data-rate  single-tone  DQPSK).  It  is  not  clear  that  the  frequency  that’s  best  for  the 
ALE  modem  will  also  be  best  for  a  Navy  data  modem.  Early  ALE  on-the-air  testing 
(DCA,  1988)  suggests  the  possibility  of  such  a  performance  disparity  between  ALE  and 
the  Harris  RF3466,  a  39-tone  DQPSK  modem. 

The  military  standard  does  not  fully  specify  the  frequency  selection  algorithm.  The 
standard  provides  the  potential  to  exchange  both  SNR  data  and  supplementary  channel 
quality  data.  The  standards  developers  recognized  the  potential  future  need  to 
exchange  channel  measurement  data  such  as  multipath  spread  and  Doppler  spread. 

OBJECTIVE 

The  objective  of  the  effort  described  was  to  identify  a  small  set  of  parameters  to 
augment  SNR  for  link  quality  assessment;  and  a  link  quality  measure  based  on 
received  SNR  plus  those  parameters,  such  that  ALE  will  select  the  best  frequency  for  a 
given  Navy  modem.  We  anticipated  that  the  link  quality  algorithm  will  be  a  function  of 
which  modem  is  to  operate  over  the  selected  frequency.  A  further  objective  is  these 
parameters  should  be  measurable  using  the  ALE  hardware  and  waveforms. 

APPROACH 

The  investigation  relied  upon  existing  analytic  results  for  the  modulations  imple¬ 
mented  in  existing  Navy  modems  and  for  the  modulation  defined  for  ALE.  Due  to 
time  and  funding  limitations,  we  further  narrowed  the  scope  to  focus  on  Navy  parallel 
tone  DQPSK  modems.  We  considered  two  channel  parameters;  root-mean-square  (rms) 
Doppler  spread  and  rms-delay  spread.  These  parameters  were  chosen  because  they 
have  a  strong  connection  to  the  analytical  results,  relauvely  simple  techniques  have 
been  developed  for  estimating  them  on  real  channels  (Perl  and  Kagan,  1986),  and  they 
can  be  easily  accommodated  within  the  ALE  data-transfer  structure  already  proposed 
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for  exchanging  channel  information.  Existing  analytic  results  were  modified  to  describe 
the  Navy  and  ALE  modems. 

This  report  presents  the  results  of  the  analytic  investigation.  Future  work  will  com¬ 
pare  the  analytic  results  with  simulation  results  obtained  using  a  Watterson-type 
HF  channel  simulator. 

The  analytic  modem  performance  models  were  examined  to  identify  benefits  avail¬ 
able  from  incorporating  rms-delay  spread  measurements  in  the  ALE  frequency  selec¬ 
tion  process. 


ANALYTICAL  ASSESSMENT 


Although  modem  performance  has  been  the  subject  of  analyses  for  at  least  four 
decades,  there  still  is  no  expression  for  performance  of  even  the  most  common  modu¬ 
lation  schemes  incorporating  a  set  of  channel  and  hardware  performance  parameters. 
However,  many  of  these  parameters  are  treated  independently.  The  functional  form  of 
the  performance  for  the  common  modulation  types  has  been  derived.  We  will  not 
attempt  to  address  performance  through  new  mathematical  analyses.  Rather,  we  will 
exploit  existing  analyses  to  build  a  plausible  model  for  empirical  performance  estima¬ 
tion. 

GAUSSIAN  NOISE 

The  simplest  analysis  concerns  binary  modem  performance  in  the  presence  of  addi¬ 
tive  stationary  white  Gaussian  noise  (Schwartz,  Bennett,  and  Stein,  1966).  There  are 
two  forms  for  bit  error  probability  (Pb)  expressions  covering  modulations  of  interest  to 
us.  The  first  of  these  is 


Pb  =  -expI-kR]  (1) 

2 

where  R  is  the  signal  energy-to-noise  ratio,  R  =  Eb/No .  Theoretical  analyses  give  k  =  1 
for  binary  DPSK  and  k  =  1/2  for  noncoherent  FSK  in  this  expression.  The  second  form 
of  the  error  probability  expression  is 


Pb  «  -^erfc[(kR)’/2]  .  (2) 

2 

Here  erfc  is  the  complementary  error  function.  The  theoretical  analyses  show  k  =  1/2 
for  coherent  FSK,  and  k  =  1  for  coherent  PSK  in  this  expression. 
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For  m-ary  DPSK,  the  symbol  error  probability  (Pe)  expression  is 


Pe  «  erfc[(kR)^/^]  . 


(3) 


An  excellent  approximation  for  m-ary  DPSK,  due  to  Arthurs  and  Dym  (1962),  can  be 
found  by  using 


k  « 


(4) 


in  equation  3. 

In  each  of  these  cases,  increasing  k  corresponds  to  shifting  the  (log-log)  error  prob¬ 
ability  curve  to  the  left  (improved  performance).  Similarly,  decreasing  k  shifts  the 
curve  to  the  right  (degraded  performance).  Such  a  shift  to  the  right  could  correspond 
to  hardware  implementation  losses  of  various  types. 

GAUSSIAN  NOISE  WITH  FLAT  FADING 

Next  consider  the  case  of  slow  flat  fading.  The  fading  is  termed  flat  because  the 
entire  signal  bandwidth  fades  together.  It  is  termed  slow  because  fading  effects  are 
assumed  constant  over  a  symbol  duration  (or  two  symbols  for  differential  phase-shift 
keying).  In  this  case,  the  following  expressions  are  well  known  (Schwartz,  et  al.,  1966). 
We  have 


Pe 


1 

2(l  +  kR) 


(5) 


for  noncoherent  FSK  and  binary  DPSK,  and 

P  «  ifi - ^  _  1 

2l  (1  +  l/kR)^/2j 

for  coherent  FSK,  coherent  PSK,  and  quartemary  DPSK.  In  these  formulas,  R  is  the 
average  received  value  of  Eb/No  ;  R  =  Eb/No  . 

Kam  (1991)  points  out  that  this  traditional  analysis  applies  to  the  case  of  zero¬ 
fading  bandwidth.  That  is,  the  fading  process  results  in  symbol-to-symbol  channel 
complex  gain  changes,  but  is  slow  enough  to  be  treated  as  constant  over  the  symbol 
duration.  For  DPSK,  this,  in  fact,  means  slow  enough  to  be  treated  as  constant  over 
symbol  pairs.  Kam  (1991)  extends  the  existing  analysis  to  2,  4,  and  8  DPSK  with 
arbitrary  order  of  diversity.  He  treats  the  case  where  the  fading  process  is  constant 


over  the  symbol  duration  but  not  across  symbol  pairs.  The  results  of  particular  interest 
to  us  are  the  case  of  no  diversity  with  binary  or  quartemary  DPSK. 

For  binary  DPSK, 

where  p  is  the  normalized  autocorrelation  function  of  the  fading  process.  For  a  Gauss¬ 
ian  fading  process,  p  is  given  by 

p  =  exp[-  (jrBT)^],  (8) 

where  B  is  the  fading  bandwidth  (b  =  for  a  Gaussian  process  with  variance  a^) 

and  T  is  the  modulation  symbol  duration.  When  p  =  1  (B  =  0)  this  expression  reduces 
to  the  expression  for  P^  of  DPSK  given  in  equation  5. 

For  quaternary  DPSK,  we  have 


This  result  applies  to  the  slow  fading  channel,  for  which  p  =  1  -  e ,  with  0  <  e 
«  1.  We  can  rewrite  this  expression  as 


We  can  now  find  the  irreducible  error  bound  as  the  limit  of  this  expression  as 
R  ->  00 .  The  result  is 

Pb  =  -j^l-- - ;; - - qjTjl.  =-{!-(!- 2£))  =  e.  (11) 

We  can  also  evaluate  the  expression  for  Pi,  in  the  case  of  zero-fading  bandwidth 
(P  =  1). 
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This  expression  can  be  compared  to  the  expression  given  in  equation  6.  In  particu¬ 
lar,  for  R  »  1/2,  equation  12  becomes 


which  corresponds  to  equation  6  with  k  =  1/4.  This  value  of  k  is  quite  close  to  the 
Arthurs  and  Dym  (1962)  value  of 

k  =  sin^l^^j  ==  0.278  .  (14) 

EXTENSIONS  FOR  DQPSK 

At  HF,  transmitted  signals  may  be  refracted  from  various  physical  locations  in  the 
ionosphere.  This  results  in  time  dispersion  of  the  received  signal.  This  time  dispersion 
is  termed  multipath.  It  produces  the  phenomenon  of  frequency  selective  fading. 

Garber  and  Pursley  (1988)  have  analyzed  the  case  of  frequency  selective  fading  for 
binary  DPSK  under  the  particular  circumstances  where  most  (>90  percent)  of  the 
delayed  energy  affects  only  the  next  symbol,  where  the  dispersion  takes  on  one  of  four 
continuous-time  distributions,  and  for  three  types  of  DPSK  waveforms.  Garber  and 
Pursley  (1988)  show  the  bit-error  performance  under  these  conditions  is  bounded  below 
by  the  performance  for  slow  flat  fading  for  low  SNR.  At  high  SNR  there  is  an  irreduc¬ 
ible  error  bound  due  to  the  intersymbol  interference.  The  principal  factors  influencing 
Pi  (the  irreducible  bound)  are  (a)  the  normalized  multipath  spread,  and  (b)  the  shape 
•  of  the  communication  signal  pulse  or  symbol.  There  is  little  performance  difference 

among  Gaussian,  exponential,  triangular,  or  rectangular  continuous-delay  distributions. 
No  analysis  has  been  identified  indicating  the  effects  of  discrete-delay  distributions  or 
of  multiple,  disjointed  continuous  distributions.  The  following  discussion  summarizes 
the  Garber  and  Pursley  (1988)  results. 

Given  the  normalized  multipath  spread  fiy,  we  can  estimate  Pi .  Bounds  for  this 
estimate  are  plotted  in  Garber  and  Pursley’s  (1988)  figure  6.  A  straight  line  approxi¬ 
mation  to  this  figure  (based  on  their  table  IV)  is 


5 


Pi  =  0.475 


(15) 


Having  found  Pj,  we  can  find  an  expression  for  Pj,.  This  expression  defines  a  curve 
that  transitions  between  the  Pb  of  the  slow  flat  fading  case  at  low  SNR  and  the  con¬ 
stant  value  Pi  at  high  SNR.  Garber  and  Pursley  (1988)  do  not  give  a  closed-form 
expression  for  this  performance  curve.  However,  the  performance  curve  is  similar  to 
that  given  by  Kam  (1991)  for  the  nonselective  fading  case.  We  will  proceed  to  con¬ 
struct  a  model  for  empirical  use  of  the  form  derived  by  Kam  (1991).  This  is  not  a 
mathematical  derivation. 


Let  g(^)  be  the  power  density  spectrum  of  the  channel  delay.  Then  the  rms  delay  is 


■[) 


and  the  rms-multipath  spread  is  defined  as 

Ai  =  NPr  ,  (17) 

where  T  is  the  communications  pulse  duration.  The  normalized-multipath  spread  Pv 
(that  is,  normalized  to  an  “effective”  communications  pulse  duration)  depends  upon 
the  shape  of  the  transmitted  communications  signal  pulses.  For  rectangular  DPSK 
pulses,  pf  =  p=  M/T.  For  sinusoid-shaped  DPSK  pulses,  Ps  can  be  found  from  the 
power  density  spectrum  of  the  channel  delay.  An  approximation  to  the  tabulated  values 
given  by  Garber  and  Pursley  (1988)  is 

Ps  =  15.1(M/T)3  °^  .  (18) 

In  total,  a  three-step  process  is  needed  to  estimate  Pi.  First,  the  channel-rms  delay 
must  be  determined.  From  this  the  appropriate  normalized-multipath  spread  must  be 
estimated,  incorporating  both  the  shape  and  duration  of  the  transmitter  symbol  modula¬ 
tion.  Finally,  the  estimated  Pi  is  calculated  using  the  normalized-multipath  spread. 

We  can  now  consider  modifications  to  the  expression  for  Pb  in  the  flat-fading  case 
to  include  multipath  frequency-selective  fading  effects.  In  particular,  multipath  results 
in  an  irreducible  error  rate  similar  to  that  observed  for  fading.  If  in  the  expression  for 
Pb  we  insert  a  term  ^  multiplying  p,  and  similarly  constrain  ^  to  the  range  0  <  ^  s  1, 
then  ^  will  have  the  same  effect  as  p.  That  is,  ^  will  introduce  an  irreducible  error 
bound  as  R  oo . 

The  resulting  expression  for  binary  DPSK  is 


Ph  =  - 


ritd-POR-i 

l  J 


Now  for  R  ->  00  and  p  =  1  (zero  fading  bandwidth  on  individual  paths)  we  have 

Pb  =  -(1  -  O  =  Pi  •  (20) 

2 

From  this  expression  we  find  the  appropriate  value  for  ^  to  be  ^  =  1  -  2Pi . 

Without  justification,  but  with  analogy  to  the  flat- fading  case,  we  model  quaternary 
DPSK  performance  as 


Let  the  multipath  have  Gaussian  distribution  with  density 

g(l)  =  exp[-^V  2a^]-  (22) 

The  “rms-multipath  delay”  is  just  the  standard  deviation  of  g(^).  The  symbol  M  is 
usually  used  to  denote  flie  rms-multipath  delay,  so  in  the  Gaussian  case  M  =  a.  The 
rms-multipath  spread  measures  the  time  spread  of  the  received  signal  relative  to  the 
symbol  duration  in  “a -units.”  These  a -units  are  multiples  of  the  standard  deviation  of 
the  multipath  spread  treated  as  a  density  function.  For  the  Gaussian  density  given 
above,  we  then  know  that:  for  =1,  about  68  percent  of  the  spread  signal  is  con¬ 
tained  within  ±1  symbol;  for  //  =  1/2,  about  95  percent  of  the  spread  is  contained  in 
±1  symbol;  and,  for  ^  =  1/3,  about  99.74  percent  of  the  spread  is  contained  in  ±i 
symbol. 

EXTENSIONS  FOR  FSK 

For  FSK  signaling,  we  can  apply  another  analysis  of  Garber  and  Pursley  (1989).  In 
general,  the  multipath  interference  results  in  an  irreducible  error  bound.  This  bound  is 
a  bit-error  probability  that  cannot  be  improved  by  increasing  the  signal-to-noise  ratio. 

At  high  SNR,  the  bit-error  rate  is  a  constant  value,  P,.  At  low  SNR,  the  bit-error  rate 
is  dominated  by  the  error-rate  curve  for  nonselective  Rayleigh  fading.  A  smooth  transi¬ 
tion  occurs  between  these  regions. 

An  estimate  for  the  value  of  Pi  can  be  obtained  using  the  Garber  and  Pursley 
(1989)  results.  This  involves  estimating  the  normalized  rms-multipath  spread 

1 

//^(h)  =  2jg'(x)x^[l  +  sinc^(hx)  -  2cos(7rhx)  sinc(hx)]dx, 

0 


(23) 


where,  for  a  power  density  spectrum  of  the  Gaussian  delay  of  the  fading  process,  we 
have 


°  (MA-hS  ^^Pl-^V2(Wr)^)  .  (24) 

In  the  expression  for  //^(h) ,  the  sine  function  is  just  sinc(y)  =  sin(;ry)/7ry.  We  see 
that  //(h)  depends  upon  the  FSK  modulation  index  h  =  2  fj  T,  where  T  is  the  baud 
duration  and  2  fd  is  the  tone  spacing,  and  upon  the  rms-multipath  spread. 

This  result  applies  to  phase-continuous  FSK.  Phase  discontinuities  result  in  irreduc¬ 
ible  error  bounds  that  are  orders  of  magnitude  higher.  Since  the  ALE  standard  modem 
is  phase  continuous,  we  will  not  consider  discontinuous  phase  FSK. 

The  integral  may  be  evaluated  numerically.  For  the  Gaussian  delay  Sj^ectrum  case, 
Garber  and  Pursley  (1989)  have  tabulated  values  of  //(h)  for  h  =  1,  2,  and  10.  This  is 
a  rather  complicated  integral  to  evaluate  since  the  integrand  is  close  to  zero  and 
involves  differences  of  values  that  are  close  to  1.  However,  we  can  approximate  the 
Garber  and  Pursley  (1989)  tabulated  data  with  the  expression 

//(h)  «  min[2.6h(MA^)^•®^  M/T]  .  (25) 

This  expression  is  plotted  in  figure  1  against  M/T,  w*th  h  as  a  parameter.  The 
tabulated  values  from  Garber  and  Pursley  (1989)  are  also  shown  on  the  figure.  In  this 
approximation,  the  exponent  may  be  changed  from  1.86  to  2  (and  the  multiplier 
increased  from  2.6  to  5.4)  with  some  decrease  in  accuracy  (errors  as  large  as  50  per¬ 
cent).  This  only  makes  sense  where  the  simpler  expression  is  required  for  ease  in 
computation. 

Given  //(h),  Garber  and  Pursley  (1989)  give  bounds  on  P,  in  their  figure  4.  We  can 
approximate  Pi  by 


Pi  *  0.29// ^h)  .  (26) 

The  functional  form  of  the  bit-error  rate  curve  for  FSK  should  be  similar  to  that 
given  by  Kam  (1991)  for  DPSK,  since  the  flat-fading  performance  form  for  FSK  is 
similar  to  that  for  DPSK.  (For  flat  fading,  the  form  for  FSK  differs  from  DPSK  by  a 
factor  of  2  dividing  the  average  SNR.)  We  then  can  model  the  bit-error  rate  perform¬ 
ance  as 


P  =  ifl  (1  -p)R/2' 
21  1  +  R/2 
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Logio  [M/n 


Figure  1.  Normalized  rms-multipath  spread  approximation. 


If  p  =  1  -  e ,  we  have 

lim  [Pb]  =  e/2  =  Pi  . 

R  -*  00 

Using  the  approximations  based  on  Garber  and  Pursley  (1989),  we  have 


e  =  0.58  min2[2.6h(M/T)’•®^M/T]. 


(28) 


(29) 


We  note,  for  small  p  (large  e),  Pb  is  close  to  1/2  for  all  R  This  does  not  corre¬ 
spond  to  a  useful  channel.  If  we  now  consider  the  crossover  point  between  the  two 
expressions  in  equation  25,  we  can  solve  to  find  the  crossover  to  be 


M/T  = 


1 

(2.6h)*-^^  ■ 


(30) 


The  largest  possible  M/T  crossover  occurs  for  h  =  2  in  the  ALE  case,  which  implies 
that  ;^(h)  =  M/T  (independent  of  h)  whenever 

NOT  >^^-0.15.  (31) 


We  can  extend  the  binary  FSK  analysis  to  8-ary  ALE  FSK.  Let  the  eight  frequen¬ 
cies  used  by  ALE  be  labeled  fi  through  fg,  as  shown  in  table  1.  Next,  assume  the 


9 


Garber  and  Parsley  (1989)  analysis  holds  for  any  pair  of  these  frequencies,  fi  and  fj. 

In  particular,  if  fj  is  transmitted,  then  the  probability  that  the  receiver  will  detect  fj  will 
follow  the  functional  form  for  binary  FSK.  However,  since  there  are  eight  tones  rather 
than  two,  the  coefficient  should  be  1/8.  Then  in  the  limiting  case  R  -»  0,  all  tones 
become  equally  likely  at  the  receiver  output,  and  P(fj  detected  |  fj  transmitted)  -»l/8. 
We  can  write 


P(fj  detected  \  fj  transmitted)  = 

Table  1.  ALE  frequencies  and  bit  representations. 


Frequency 

Hz 

Bit  Pattern 

fi 

750 

000 

1000 

001 

fs 

1250 

on 

f4 

1500 

010 

fg 

1750 

no 

^6 

2000 

111 

f7 

2250 

101 

fs 

2500 

100 

Here  we  assume  that  R  =  Eb/No  is  the  average  received  energy  per  bit  divided 
by  the  noise-spectral  density,  and  3R  is  the  average  received  energy  per  symbol  (tri-bit) 
divided  by  the  noise-spectral  density.  We  also  have 


P  =  1  -  0.58  min2[2.6h(M/T)^  ®MVI/T] 


(33) 


and  h  =  2fd  T  =  1  fj  -  fj  1  T.  Since  for  ALE  we  have  T  =  1/125  second,  we  can  write 
h  =  2|i-j|. 

The  number  of  bit  errors  given  an  incorrect  frequency  decision  at  the  receiver 
depends  on  the  specific  values  of  i  and  j.  Table  2  shows  the  number  of  bit  errors  as  a 
function  of  ji  -  j|.  If  we  assume  all  bit  patterns  are  equally  likely  at  the  source,  we  can 
write  the  received  bit-error  probability  as 


Pb 


SS 


P(fj  detected  |  fj  transmitted) 


i  =  1  j  s-s  i 

Here  Njj  is  the  number  of  bits  of  the  tri-bit  in  error  when  fj  is  detected,  given  that  fj  is 
transmitted. 
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Table  2.  Possible  error  occurrences. 


1  i-j  1 

bit  errors 

1 

1 

2 

2 

3 

3  single,  2  triple 

4 

2 

5 

1  single,  2  triple 

6 

2 

7 

1 

For  small  values  of  p  (£  close  to  1),  meaning  for  large  M/T,  p  is  independent  of  h 
and  hence  independent  of  i  and  j.  This  independence  occurs  for  M/T  satisfying 
(M/T)°  >  1/(2. 6h).  For  ALE,  p  is  independent  of  i  and  j  whenever  M/T  >  0.15, 

since  the  minimum  possible  h  is  2.  This  occurs  whenever  p  <  0.987.  In  this  case  we 
have 


Pb 


P(fj  detected  |  fj  transmitted) 
8 


S 

i  =  1 


I 

j  i 


_L  r  1  +  (1  -  P)3R/2] 

64  L  1  +  3R/2  J  Zj  ZL  3 

i  =  1  j  ^  i 

ri  +  (1  -  P)3R/2l 
2  L  1  +  3R/2  J 

1  ri  +  0.58(M/T)23R/2‘1 

2  L  1  +  3R/2  J 


In  the  case  of  R-»  oo,  the  irreducible  error  bound  is  given  by 


(34) 


P,  =  0.29(Mn')2  , 


(35) 


which  applies  whenever  M/T  >  0.15.  The  value  of  Pj  for  M/T  =  0.15  is  Pj  =  0.0065- 

We  can  also  find  the  bit-error  probability  expression  for  large  values  of  p  (e  close 
to  zero).  This  corresponds  to  (M/T)°  ®®  <  1/(2. 6h).  For  ALE  the  maximum  value  of  h 
is  14,  so  “large  p”  corresponds  to  M/T  <  0.015.  In  this  case  we  have 
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8 


Ph  = 


J_  'V'  ^  [1  +  0  -  P)3R/2'j  ^ 
64  2-ui  2^  L  1  +  3R/2  J  3 


i  =  1  j  ^  i 

1 


64(1  +  3R/2)  2  ^  ^ 

i=  1  j  !“  i 
1^32  +  7.84R(M/r)3-^2 


(36) 


1 


64(1  +  3R/2) 


i  =  1  j  5^  i 


1 


^  [1  +  310R(MAr)3-’2  3  _ 


2(1  +  3R/2) 

Here  we  have  explicitly  shown  the  dependence  of  £  on  h,  and  hence  on  i  and  j,  by 
writing  ey .  As  R  -*  * ,  the  corresponding  irreducible  error  bound  is 


P,  =  103(M/T) 


3.72 


(37) 


The  value  of  Pi  for  M/T  =  0.015  is  Pi  =  1.7  x  10*^. 

The  use  of  the  “large  e”  approximation  provides  an  upper  bound  to  Pi  whenever  e 
is  not  “large”  (smaller  M/T).  In  particular,  when  M/T  =  0.015  the  “large  e”  approxi¬ 
mation  gives  Pi  =  6.5  x  10'^,  about  four  times  the  value  determined  from  the  “small 
f  ”  approximation.  Since  even  this  larger  error  rate  is  quite  small  for  practical  HF 
systems,  we  may  safely  use  the  “large  e”  approximation  in  all  cases. 


ANALYTIC  ASSESSMENT  CONCLUSIONS 


Suppose  an  HF  channel  has  delay  spread  but  no  Doppler  spread.  We  can  then  use 
the  FSK  and  DQPSK  expressions  and  compare  modem  performance  for  this  channel. 
Figures  2  and  3  show  the  bit-error  performance  for  DQPSK  and  FSK  respectively.  We 
assumed  2400  bit/second  data  rate  via  a  16-tone  modem  with  no  error  correction  cod¬ 
ing  for  the  DQPSK  case  of  figure  2.  In  figure  3  we  used  the  parameters  of  an  ALE 
modem  in  estimating  FSK  performance.  The  analysis  leading  to  figure  2  applies  only 
for  multipath  confined  to  one  symbol,  or  M  <  125  (log  M  <  2.09).  Figure  3,  therefore, 
shows  an  error  (Pi,  >  1/2)  for  large  M.  The  true  bit-error  probability  approaches  (does 
not  exceed)  1/2  for  large  M. 
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The  first  point  to  note  about  these  figures  is  that  they  are  very  similar.  This  simi¬ 
larity  implies  that  performance  of  the  ALE  modem  should  be  a  useful  guide  for  pro¬ 
jecting  DQPSK  modem  performance  on  this  type  of  channel. 

The  second  point  to  note  is  that  modem  performance  is  strongly  dependent  upon 
the  rms-multipath  spread,  M.  SNR  alone  is  not  a  reliable  predictor  of  bit-error  per¬ 
formance.  However,  since  SNR  and  Pb  together  may  be  used  to  estimate  M  for  this 
channel,  SNR  could  be  used  with  either  Pb  or  M  to  select  a  best  frequency.  In  practice, 
this  is  not  the  case  since  an  accurate  estimate  of  Pb  is  difficult  to  obtain,  especially 
since  only  a  relatively  few  bits  are  exchanged  in  the  ALE  algorithm.  Furthermore, 
Doppler  spread  may  also  significantly  affect  modem  performance— then  all  three 
parameters  (SNR,  Doppler  spread,  and  delay  spread)  are  needed,  and  SNR  plus  bit¬ 
error  rate  are  insufficient  for  predicting  performance. 

The  DQPSK  analysis  shows  Doppler  spread  and  delay  spread  are  both  significant 
parameters  for  predicting  modem  performance.  The  implication  of  the  analysis  is  a 
technique  for  measuring  delay  and  Doppler  spreads  should  be  incorporated  into  ALE 
handshaking,  and  these  values  together  with  SNR  should  be  used  in  the  frequency- 
selection  algorithm.  The  frequency-selection  algorithm  would  then  take  modem  charac¬ 
teristics  into  account  when  prioritizing  the  available  frequencies. 


14 


REFERENCES 


Arthurs,  E.,  and  H.  Dym.  December  1962.  “On  the  Optimum  Detection  of  Digital  Sig¬ 
nals  in  the  Presence  of  White  Gaussian  Noise,  TRE  Trans.  Comm.  Sys.,  Vol.  CS-10, 
No.  4,  pp.  336-372. 

DCA  Joint  Tactical  Command,  Control,  and  Communications  Agency,  88-DCA-T002. 
August  1988.  HF  ALE  Equipment  Performance  Test  Report,  by  Gail  Query,  pp.  B-88 
to  B-89. 

Garber,  F.  D.  and  M.  B.  Pursley.  January  1988.  “Performance  of  Differentially  Coher¬ 
ent  Digital  Communications  Over  Frequency  Selective  Fading  Channels,”  IEEE 
Trans,  on  Comm.,  Vol.  36,  No.  1,  pp.  21-31. 

Garber,  F.  D.  and  M.  B.  Pursley.  January  1989.  “Performance  of  Binary  FSK  Commu¬ 
nications  Over  Frequency-Selective  Rayleigh  Fading  Channels,”  IEEE  Trans,  on 
Comm.,  Vol.  37,  No.  1,  pp.  83-89. 

Kam,  P.  Y.  February  1991.  “Bit  Error  Probabilities  of  MDPSK  Over  the  Nonselective 
Rayleigh  Fading  Channel  with  Diversity  Reception,”  IEEE  Trans.  Comm.,  Vol.  39, 
No.  2,  pp.  220-224. 

Perl.  J.  M.  and  D.  Kagan.  January  1986.  “Real-Time  HF  Channel  Parameter  Estima¬ 
tion,”  IEEE  Trans.  Comm.,  Vol.  COM- 34,  No.  1,  pp.  54-58. 

Schwartz.  M.,  W.  R.  Bennett,  and  S.  Stein.  1966.  Communication  Systems  and  Tech¬ 
niques,  p.  406,  McGraw-Hill. 


15 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  Information  Is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  Instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  Informatt^  Operations  and  Reports.  1 21 5  Jefferson  Davis  Highway.  Suite  1 204.  Arlington.  VA  22202-4302. 
and  to  the  Office  of  Management  and  Budget  Paperwork  Reduction  Pro|eci  (0704-0188)  Washington.  DC  20503 


2  REPORT  DATE  3  REPORT  TYPE  AND  DATES  COVERED 


1  AGENCY  USE  ONLY  (Leave  blank) 


December  1991 


4  TITLE  AND  SUBTITLE 


ALE  FREQUENCY  SELECTION  FOR  NAVY  MODEMS:  ANALYSIS 


6  AUTHORIS) 

C.  Fuzak,  Jr. 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152-5000 


5  FUNDING  NUMBERS 

PR:  ZF09 
WU:  DN300172 
PE:  0602936N 


8  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Chief  of  Naval  Research 

Independent  Exploratory  Development  Program  (lED) 

OCNR-20T 

Arlington,  VA  22217 


10  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


NOSC  TR  1461 


12a  DISTRIBUnON/AVAILABILfTY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


13  ABSTRACT  (Mmmum  200  f/norcts) 


The  objective  is  to  develop  a  simple  model  for  HF  (2-30  MHz)  modem  performance  as  a  function  of  various  channel 
parameters,  and  to  develop  techniques  for  predicting  those  parameters. 

Modem  performance  was  modeled  based  on  three  parameters— signal-to-noise  ratio  (SNR),  delay  spread,  and 
Doppler  spread.  Each  of  these  parameters  is  measurable  using  known  techniques.  The  model  demonstrates  the  inadequacy  of 
using  only  SNR  and  bit  error-rate  estimates  for  frequency  selection  in  automatic  link  establishment  (ALE)  algorithms. 

The  benefits  from  incorporating  delay  and  Doppler  spreads  in  the  ALE  frequency  selection  algorithm  are 
sufficiently  clear  to  justify  extension  of  the  algorithms  to  include  these  factors.  Further  research  is  required  to  verify  the 
sufficiency  of  SNR,  delay  spread,  and  Doppler  spread  for  frequency  selection. 


17  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIHED 


18  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19  SECURITY  CLASSIFICATION 
OF  abstract 

UNCLASSIFIED 


15  NUMBER  OF  PAGES 


16  PRICE  CODE 


20  UMITATION  OF  ABSTRACT 


SAME  AS  REPORT 


NSN  754001 -280-55C0 


Sia.'Kjard  fom  298 


NSN  754CM)1-260'5500 


Standard  form  296 

UNCLASSmED 


INITIAL  DISTRIBUTION 


Code 

0012 

Patent  Counsel 

(1) 

Code 

0142 

K.  Campbell 

(1) 

Code 

0144 

R.  November 

(1) 

Code 

761 

J.  Loughlin 

(1) 

Code 

808 

G.  Clapp 

(1) 

Code 

808 

R.  Peterson 

(1) 

Code 

824 

C.  Deneris 

(1) 

Code 

824 

P.  Donich 

(1) 

Code 

824 

C.  Fuzak 

(10) 

Code 

824 

K.  Owens 

(1) 

Code 

824 

J.  Rahilly 

(1) 

Code 

951 

T.  Danielson 

(1) 

Code 

952B 

J.  Puleo 

(1) 

Code 

961 

Archive/Stock 

(6) 

Code 

964B 

Library 

(3) 

Defense  Technical  Information  Center 
Alexandria,  VA  22304-6145  (4) 

NCCOSC  Washington  Liaison  Office 
Washington,  DC  20363-5100 

Center  for  Naval  Analyses 
Alexandria,  VA  22302-0268 

Navy  Acquisition,  Research  &  Development 
Information  Center  (NARDIC) 

Alexandria,  VA  22333 

Navy  Acquisition,  Research  &  Development 
Information  Center  (NARDIC) 

Pasadena,  CA  91106-3955 

Naval  Research  Laboratory 
Washington,  DC  20375-5000 

Space  &  Naval  Warfare  Systems  Command 
Washington,  DC  20363-5100  (3) 


1 


